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Abstract 

We consider solutions of eleven-dimensional supergravity constructed in [1, 2] that 
are half-BPS, locally asymptotic to AdS-/ x and are the holographic dual of heavy 
Wilson surfaces in the six-dimensional (2,0) theory. Using these bubbling solutions 
we calculate the holographic entanglement entropy for a spherical entangling surface 
in the presence of a planar Wilson surface. In addition, we calculate the holographic 
stress tensor and, by evaluating the on-shell supergravity action, the expectation value 
of the Wilson surface operator. 
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1 Introduction 


The celebrated proposal of Ryu and Takyanagi [3, 4] provides a method to calculate en¬ 
tanglement entropies for held theories that have holographic duals. Originally the proposal 
was used to calculate the entanglement entropy for theories in their vacuum state, but was 
quickly generalized to include more general settings, such as hnite temperature and time- 
dependent states (see e.g. [5] for a review). For spherical entangling surfaces it was observed 
by Casini, Huerta and Myers [6] that the holographic entanglement entropy can be mapped 
to the thermal entropy of a hyperbolic black hole. In the held theory, the corresponding 
entanglement entropy is mapped to the thermal entropy on a hyperbolic space. 

Another generalization concerns entanglement entropy in the presence of extended de¬ 
fects, such as Wilson lines. In the probe approximation these defects are described by strings 
or branes in the AdS space and the backreaction of the stress energy of the branes on the 
bulk geometry is neglected. For example, the holographic description of a Wilson loop in 
SU{N) Af = 4 SYM in the fundamental representation is given by a fundamental string in 
AdS^ X [7, 8], whereas higher dimensional representations can be described by D3- (D5-) 
branes with AdS 2 x worldvolume in AdS^ x [9, 10]. See [11, 12, 13] for a discus¬ 

sion of entanglement entropy in the presence of probe brane defects. When the dimension of 
the representation increases and becomes of order the backreaction cannot be neglected 
and the probe is replaced by a new bubbling geometry with flux. The bubbling holographic 
solutions corresponding to half-BPS Wilson loops in A/" = 4 SYM were found in [14, 15, 16]. 

In [17], Lewkowycz and Maldacena applied the Casini, Huerta and Myers mapping to 
the calculation of the entanglement entropy in the presence of Wilson loops in A7 = 4 
SYM theory and ABJM theories [18]. They showed that the entanglement entropy can be 
calculated from the expectation value of the Wilson loop operator as well as the one point 
function of the stress tensor in the presence of the Wilson loop. For BPS Wilson loops 
these quantities can be evaluated using localization and reduced to matrix models [19, 20]. 
In [21] by two of the present authors, the holographic entanglement entropy was calculated 
using the bubbling solution dual to half-BPS Wilson loops. It was shown that the result 
agrees with [17] when the exact map between matrix model quantities and the supergravity 
solution, found in [22, 23], is applied. 

The goal of the present paper is to generalize the holographic calculation of the entan¬ 
glement entropy to the case of six-dimensional (2, 0) theory with half-BPS Wilson surfaces 
present. The (2, 0) theory can be dehned either by a low-energy limit of type HB string 
theory on an A^^i singularity [24] or by a decoupling limit of N coincident M5 branes [25]. 
While there exists no simple Lagrangian description of this theory due to the presence of 
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Figure 1: The spherical entangling surface dA is the boundary of a region ^ on a constant 
time slice of the (2, 0) theory on M®. The Wilson surface T intersects this surface twice. 

tensor helds B~^ with self-dual held strength, the holographic dual [26] is given by M-theory 
on AdSj X S^. In analogy with the Wilson loop one expects that the six-dimensional theory 
has extended Wilson surface operators [27] of the form 

Wr ~ trexp 5+^ (1.1) 

In the probe approximation the Wilson surface operators can be described by embedding 
M2-branes [28, 29] or M5-branes [30, 31, 32] on various submanifolds inside the AdS^ x S^. 
It is an interesting open question whether the expectation value of Wilson surface operators 
can be calculated by localization in the (2, 0) theory. 

In [1, 2], bubbling solutions corresponding to half-BPS Wilson surfaces were found (see 
[33] for earlier work in this direction). In the present paper we will use these solutions 
of eleven-dimensional supergravity to calculate the entanglement entropy as well as other 
holographic observables. The bubbling solutions feature an AdSs x x fibration over 
a two-dimensional Riemann surface S with boundary. The solutions are locally asymptotic 
to AdSj X S‘^ and the six-dimensional asymptotic metric on the AdS^ boundary is naturally 
AdSs xS^. It is convenient to describe the Wilson surface on this space by imposing boundary 
conditions at the boundary of AdS^ and choosing Poincare coordinates for the AdS^ factor 
describes a planar Wilson surface. However, as discussed in section 3.3, this metric on 
AdSs X can be related to the more familiar flat metric on by a conformal transformation. 
It is easier to visualize the geometry of our setup on M®: the entangling surface at constant 
time is a four-sphere of radius R and the Wilson surface is a line (also Filing out the time 
direction) that intersects the four-sphere at two points, as illustrated in figure 1. 
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1.1 Summary of results 


For the convenience of the reader we will present our main results here and put them into 
context. The change of the entanglement entropy due to the presence of the Wilson surface 
is given by 


AS^ = 


4iV3 

IT 


16 + 3ml — Sms 
320 


i<j 



( 1 . 2 ) 


We also calculate two other holographic observables. The stress tensor in the presence of 
the Wilson surface for the AdS^ x boundary coordinates is given by 

A{Tij) dx^dx^ = (16 + 3m^ - Sms) [ds\^s^ - ds%) (1.3) 

Finally we calculate the expectation value of the Wilson surface operator by evaluating the 
regularized on-shell supergravity action: 

log(Wr) = -^Yo\{AdS,) {F + U) (1.4) 

iO/TT 

The m 2 ,s are quantities that depend on the parameters of a general bubbling solution. The 
cut-off T] is the distance from the Wilson surface, as discussed in more detail in section 3.3. 
Also, is a hnite one-dimensional integral that is dehned in section 5.3.2. 


1.2 Structure of the paper 

The structure of the paper is as follows. In section 2 we review the bubbling half-BPS solu¬ 
tions of M-theory originally obtained in [1, 2] and work out the behavior of the solution near 
the asymptotic boundary. In particular, we determine the Fefferman-Graham map for an 
asymptotic AdS^ x boundary metric. In section 3 we calculate the entanglement entropy 
for a spherical entangling surface following the Ryu-Takayanagi prescription for the bub¬ 
bling solution. In section 4 we use the methods of Kaluza-KIein holography and holographic 
renormalization to calculate the one point function of the stress tensor for the bubbling 
solution. In section 5 we evaluate the on-shell action of eleven-dimensional supergravity to 
determine the expectation value of the Wilson surface. We show that the bulk part of the 
action is given by a total derivative and evaluate the integral as well as the Gibbons-Hawking 
term. In section 6 we provide a brief discussion of our results as well as possible avenues 
for future research. In the interest of readability we present technical matters and details of 
calculations in several appendices. 
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2 Review of bubbling M-theory solutions 


In this section we will review the construction of half-BPS M-theory solutions found in [1] 
that are locally asymptotic to AdSj x S^. These solutions generalize the construction of Janus 
solutions [34, 35] in type IIB to M-theory. They correspond to the holographic description 
of Wilson surface defects in the six-dimensional (2, 0) theory, where the Wilson surface is 
‘heavy’ and the backreaction on the geometry is taken into account. 

One demands that these solutions preserve an OS'p(4*|2)©OS'p(4*|2) sub-superalgebra of 
the OS'p(8*|4) superalgebra of the AdSj x vacuum. This form of the preserved superalge¬ 
bra is uniquely determined by demanding that the solution has sixteen unbroken supersym¬ 
metries and preserves so{2, 2|M) associated with conformal symmetry on the worldvolume of 
the Wilson surface, so(4|]R) corresponding to rotational symmetry in the space transverse to 
the Wilson surface and an unbroken so(4|M) R-symmetry [36]. Note that a generalization 
was recently analyzed in [37] in which the preserved superalgebra is iJ(2|l,7) © iJ(2|l,7), 
but we will not discuss this case here. 

It follows from these superalgebra considerations that the bubbling BPS solution has an 
so(2,2|M) © so(4|]R) © so(4|]R) algebra of isometries. Furthermore, the solution preserves 
sixteen of the thirty-two supersymmetries. The BPS equations were solved in [1] and the 
global regular solutions were found in [2]. The ansatz for the eleven-dimensional metric is 
given by an AdS^ x x fibration over a Riemann surface S with boundary: 


= /f + fi + /I + V Mii| 


( 2 . 1 ) 


where we denote the complex coordinate of the two-dimensional Riemann surface by v. In 
addition, dsla and dslg are the metrics on the unit-radius three-spheres and the metric on 
the unit-radius Euclidean AdSs in Poincare half-plane coordinates is given by 

dz"^ + dt'^ + 


dA 

^^AdS3 


( 2 . 2 ) 


The Wilson surface on the boundary AdS'i x fills the t, I directions and is located at z = 0. 

These solutions are parametrized by a harmonic function h and a complex function G{v, v) 
that satisfies a first order differential equation: 


d^G= - {G + G)dAnh 


(2.3) 


It is useful to introduce the combinations^ 


W+= \G - G\ + 2\G\\ W_ = \G - G\ - 2\G\^ (2.4) 

^Note that there is a typo in eq (2.5) of [2]: it should read W'^ = —— (G — G)^. 
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Figure 2: AdS^ x parameterized on the half strip. 


in terms of which the metric functions in (2.1) are given by 


/? = ih\l - |G|2) 

, _ hHV^W. „ ^ 

16(1-|G'|2)2’ ^ 


f! 


4^2(1- |G|2) 


hF_ 

wl 


(dyhdyh)^ 

w? 


(i-i^r) 


W+W_ 


(2.5) 


It was shown in [2] that for a solution to be regular the functions h and G must satisfy the 
following conditions on the Riemann surface S and its boundary: 


h = 0, G = 0, +i, V G cIS 

/i>0, \G\^<1 ugS (2.6) 

First we consider the simplest example: the AdSj x vacuum solution. This can be obtained 
by choosing S to be the half strip = {v = p + iq/2,p > 0, g G [0, vr]} with 

h = —iL^ (cosh(2n) — cosh(2h)), G = —i -— . ^ (2.7) 

olllll ^ (J 

which produces 

fi = 2Lcoshp, /2 = 2Lsinhp, = Lsinq, p = L (2.8) 

Hence the AdS^ x metric is given by 

= 4L^ (dp^ + cosh^pds^^^g + sinh^pdsla) + (dg^ + sin^ gdsig) (2.9) 

This geometry is represented in figure 2. More general bubbling solutions can be constructed 
once we realize that the AdSj x solution can be mapped from the half strip to the upper 
half-plane via 


w = cosh(2n) 


( 2 . 10 ) 


7 













Figure 3: (a) AcISy x on the upper half-plane, (b) General bubbling solution with n 
four-cycles C^\ i = 1,... ,n and n — 1 four-cycles C^\ i = 1,... ,n — 1. 


The functions h and G then take the form 

'3/ _\ ^ ''/ tC-|-l 


h = —iL'^{w — w), 


«=2 


w — 1 


( 2 , 11 ) 


<y(w + 1)(mi + 1) \/(ui - l)(m - 1) 

Note that the boundary of S is now located at the real line and that on the real line the 
function G = +i when Rew E [—1,1] and G = 0 when Retc > 1 or Reta < —1. 

A general bubbling solution is constructed by choosing a simple form for h and the 
following linear superposition for G\ 

2n 

h=-iL^{w-w), G = g{^) =-- \ (2.12) 




i=l 


where w is now a general coordinate on the upper half-plane. The solution is completely 
characterized by the choice of 2n real numbers with i = 1, 2,..., 2n that are ordered 


-CX) = ^0 < Cl < 6 < ■ ■ ■ < 6n < C2n+1 = +00 


(2.13) 


G| 


Im iR=0 


(2.14) 


We have introduced Co and C 2 n+i to simplify the expression for the boundary condition that 
the function G satisfies on the real line: 

0 Rew e [C 2 fc,C 2 fc+i] k = 0,1,2,... ,n 

+i Rew e [C 2 fc+i,C 2 fc+ 2 ] ’ A: = 0,1,2,... ,n - 1 

A general bubbling solution is characterized by the appearance of new nontrivial four- 
cycles (see figure 3). The n four-cycles are constructed by connecting two boundary 
points on different intervals where the volume of the three-sphere shrinks to zero. This 
generalizes the construction of the four-sphere in the AdSj x vacuum solution. In addition, 
the geometry also has n — 1 four-cycles that are constructed by connecting points on 
different intervals where the three-sphere 5^ shrinks to zero size. In the bubbling solution 
these cycles carry nontrivial four-form flux and are the remnants of M5-branes wrapping 
AdS^, X S'^ and AdS^, x S'^, respectively. 










2.1 Asymptotic behaviour and regularization 


In this section we study the asymptotic behaviour of a general bubbling solution. We will 
see later that the area integral and the action integral both diverge, so we need to regulate 
the integrals and map the regulator to the Fefferman-Graham (FG) UV cut-off. 

It is convenient to choose the following coordinates on S: w = re*®. The boundary of 
AdS^ X is located at r —)■ oo. The expressions for G and G given in (2.12) can be expanded 
at large r in terms of the generating function of the Legendre polynomials 


with the result 


1 

\/l — 


fc =0 


(2.15) 


G 


i^ak{9)mk ^ 

k=\ 


i ak{ff)rak 

k=l 


where the dependence on the angular coordinate 6 is given by 


(2.16) 


ak{9) = Pk-i{cos9) - e*®Pfc(cos 6 ') 

dk{9) = Pk-i{cos9) - e"*®Pfc(cos 6 ') (2.17) 


The moments ruk are dehned via 


2n 


rrik = 




(2.18) 


i=l 


To ensure that a general bubbling solution is asymptotic to AdS^ x with radii 


(2.19) 


we must identify mi = 2 . This provides a constraint on the Also, for the AdSj x 
solution (n = 1 ) we note here that all even moments vanish and all odd moments equal 2 . 
The metric functions take the following forms as power series in large r: 

ff ^ ^ 4 — m 2 cos 9 ^ 3(8 + ml — 2m^) -|- (8 -|- Sm^ — lOms) cos 29 , ^ / 1 

Wo — 2v “t“ ~ “t“ ; “t” G I ~ 


L2 

/I . 

^ = 2r 


f: 


L2 sin^ 9 


= 1 + 


2 ' 24r 

4 -I- m 2 cos 6 ^ ^ 3(8 -I- m^ — 2 m 3 ) -I-(8-I- 3m| — lOms) cos26' 

2 ^ 247 

m 2 cos 9 3(m| — 8 m 3 ) -|- (32 -|- 3 m 2 — 40m3) cos 29 


+ 0 - 


2r 


96 r2 


po ¬ 


rn 


P_^ _ 

12 4^2 


1 m 2 COS 9 —3(16-|-ml — 8 m 3 )-|-(16 — 3m2-|-40m3) cos 20 

P ^ O P . A P O I 


8H 


384 




( 2 . 20 ) 
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Next we present the mapping of the (r, 6) coordinates for large values of r to an FG coordinate 
system (m, 9) for a general bubbling solution. We need this map to dehne the large r cut-off 
function as well as to perform the Kaluza-Klein (KK) reduction in the calculation of the 
expectation value of the stress tensor in section 4. 

It is natural to consider a Wilson surface living on AdS^ x S^. This space is related to 
by a Weyl rescaling. We can choose to adapt our FG chart to either space; here we choose 
the former. The general FG metric that preserves the AdS^ x x isometry of a bubbling 
solution is given by 


ds^ = 


—^ (^du (y.2dsg3^ -|- cx^dO -|- C(4.ds 


u 


( 2 . 21 ) 


Equating this metric with the bubbling metric (2.1) we hnd 


fl = 


AL^ai 




/! = 


4:L^a2 




fl = 


at2j 2 

4p2 (^^2 ^ ^ + L‘^a,de^ (2.22) 

We regulate the spacetime at a small value of u and identify this with e: the (dimensionless) 
UV cut-off on AdS^ x S^. The boundary conditions on the coordinate map and the functions 
«*(«, 6) at small u must be chosen to ensure that the boundary metric is ds\^g^ + ds^g^ and 
the transverse is recovered. We hnd 


r = -^-(-..., 0 = 9 + ... 

oi\ = 1 -|- ..., 0,2 = 1 “ 1 “ ..., 03 = 1 -|- ..., 04 = sin 9 . (2.23) 


Whilst we have not been able to solve (2.22) in closed form, we can build the coordinate 
map as an asymptotic expansion in u. The mapping is given by 


and 


r = ^ + 


1712 COS 9 3(—16 — 1712 + 8 ^ 3 ) + (16 — 21 m 2 -l- dOrus) cos 29 


+ 




+ 


768 


cos^ 

18432 


(48m2 — 43m2 -l- 40m2m3 -|- 80m4 


— (48m2 — 203m2 -|- 680 m 2 m 3 — 560m4) cos26'j -|- 0{ 




~ m 2 sin 0 2 (16 — 27m2-l-40m3) cos6'sin6' . 

9 = 9 --- u -——- u 


768 


sin^ 

18432 


(296m2m3 — 48m2 — 98m2 — 200m4 


(48m2 — 139m2 -|- 400m2m3 — 280 m 4 ) cos 20 j u® -|- 0 (m®) 


u 


(2.24) 


(2.25) 
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The area integral and action integral both diverge at large r. It is useful to express the 
coordinate map as a cut-off relation Vc = rc{0, e). This is found by hrst inverting the relation 
(2.25) in the small u limit and then eliminating 6 from (2.24). The result is 




1712 COS 9 —3(16 -I- 5m2 — Sm^) + (16 — 9ml + dOms) cos 29 


+ 


768 


+ 


cos 9 
9^ 


(—48m2 -|- 557712 — 16077727773 -|- 407774 


(487772 -I- 257772 — 16077727773 -|- 280?774) COS 29) -|- 0{e^) 


(2.26) 


3 Holographic entanglement entropy 


The Ryu-Takayanagi prescription [3, 4] states that the entanglement entropy of a spatial 
region A is given by the area of a co-dimension two minimal surface At in the bulk that is 
anchored on the AdS boundary at dA: 

A . 

-^min 


^.4 = 


AG 


( 11 ) 

N 


(3.1) 


Since we are dealing with static states of our CFT, this surface lies on a constant time slice. 
If this surface is not unique, we choose the one whose area is minimal among all such surfaces 
homologous to A.^ 

In the following section we derive the minimal surface At for a general bubbling solution 
and show that its restriction to the boundary maps to a four-sphere in M®. We then evaluate 
its regulated area and compare with our expectations from M®. 

3.1 Minimal surface geometry 

A bubbling geometry is an AdS^ x x hbration over E. We consider a surface At at 
constant t that hlls the x and has prohle 2 ; = z{w,w,l), where is the AdS^ radial 
coordinate dehned in (2.2). The area functional becomes 


dz 


The equations of motion derived from this functional are solved by 

z{w^w^lY + 1“^ = (3.3) 

This semicircle is simply a co-dimension two minimal surface in AdS^. Following [12, 13] it 
is straightforward to see that this is in fact the surface of minimal area (within this ansatz). 


^This minimal surface prescription was recently established on a firm footing by the analysis of 


11 










The surface (3.3) is independent of the AdS^ radial coordinate. Thus, the boundary dA 
of the entangling region on AdSs x is given by the same formula. To understand this, let 
us consider two coordinate charts on M®: 

ds^6 = = dt^ + [dx^ + sin^ y ds^s) (3.4) 

The map between these two charts is given by 

z = fsmx, l = fcosx (3.5) 

Thus, our dA on AdS^ x can be written as a four-sphere of radius R on (given by 
f = R) after a Weyl rescaling. 


3.2 Evaluating the area integral 

The combination of metric factors that appears in the area integral (3.2) can be written 


hflfip" = \\d^h\^hW. 

The entanglement entropy is proportional to the area evaluated on the surface (3.3): 

^ VoliS^Y R 


(3.6) 


where we have dehned 


8G 


( 11 ) 

N 


R^-P 


(3.7) 


Ji = / d‘‘w\dji\^ h |G - G| 

Jt, 

J 2 = -2 / ifw \a„,h\^ h \Gf‘ 


(3.8) 


Substituting w = r e™ into (2.12) we hnd for Jj 


1 


Ji =-4L9 V(-1)M dOsmO 
Jo Jo 


r^(r cos 9 — ^i) 


rc(9,£) 

dO sind / dr , _ 

0 io ,/r7 + - 2r^i cos 0 


(3.9) 


The overall minus sign follows from the fact that G — G < 0 on the upper half-plane. We 
carefully evaluate this expression in appendix A.l. The hnal result is given in equation (A. 5) 
and takes the form 


Ji = L® 


64 —24 -|- 3m2 — Sms 


3^4 


+ 


15 


+ o (e‘) 


(3.10) 
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Next we consider the second term 


J 2 = —2 / (fw \dujh\‘^ h |Gp 

is 

r-K rrc{9,e) 

= -2L^ / de sine 


dr 


X 


2n + 2^(-l) 


i+j _ 


-rcos9{^i + Q + 


i<j \/r‘^ - 2r^i cos 0 + - 2r^j cos 6* + 

We carefully evaluate this integral in appendix A.2 and the hnal result is (A.22) 

64 4 


(3.11) 


J2 = L^ 


3e2 3 


5 ^(-i)*+^'ie«-e,f+ 0(5') 


i<j 


(3.12) 


Note that the second term cannot be expressed in terms of the moments ruk- 

Now we handle the integral over 1. Recall that the minimal surface formula (3.3) describes 
a semicircle for which z G [0, R] and I G [—R, R]- Note that Ji ,2 are independent of I because 
the cut-off function is. The I integral diverges at both limits; rewriting via (3.3) as an integral 
over z, we regulate with a cut-off at x = r/: 

D rR 


dl 

^-\/R^-v 


2 R^-P 


= 2 


= 21og 


dl 


R 


R^-P 

’R+^/R^-t]^ 
V 


= 2 


dz 


R 




= 2 log — 


(?)- 


rj 


2R2 


+ 0 ( 77 ") (3.13) 


Finally we put these pieces together to compute the divergent entanglement entropy 
(3.7): 


R. = 


L^Vol{S' 


3\2 


4G 


(11) 

N 


64 64 ^ —24 -|- 37772 — 87773 

3 ? ~ 3 ? 15 




i<j 




Employing the dehnitions 

L = (tt ip, = 2 Vol(5^) = 277^ 

this becomes 


(3.14) 


(3.15) 


5.1 = 


4A1S 


1 1 —24 -|- 37770 — 87773 


£4 £2 + 


320 




i<j 




(3.16) 
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Evaluating this result on the vacuum we hud 


qiO) 


4iV3 


' 1 1 3 ^ , 2 ^ 


log(— 
V V 


(3.17) 


Subtracting this vacuum contribution from (3.16) we arrive at our hnal result for the change 
in entanglement entropy due to the presence of the Wilson surface: 




4iV^ 16 + 3m2 — 8m3 
IT 320 


i<j 



(3.18) 


Note that the power divergences with respect to the EG cut-off e are cancelled in this 
subtraction and only a logarithmic divergence in r] remains. 


3.3 Physical interpretation 

In this section we give a physical interpretation of our result for the entanglement entropy. 
First, recall that during the calculation we introduced two separate regulators. The EG cut¬ 
off £ can be viewed as a regular UV cut-off for a holographic theory with a six-dimensional 
AdSs X boundary. In addition, when performing the integral over the AdS^ coordinate I 
in (3.13) we introduced a cut-off r] on the AdS^ radial coordinate z in Poincare slicing (2.2). 
One might be tempted to view r] as purely an IR cut-off that regulates the inhnite volume of 
the boundary theory. However, it also has an interpretation as a UV cut-off on the minimal 
distance to the Wilson surface in the boundary theory. 

This interpretation is most easily demonstrated by considering the vacuum spacetime. 
The map that relates AdS^ x with AdS^ x boundary in (2.9) to a metric with 
boundary 

4L^ 

{dv? + dt^ + df -I- -|- r'^ds'^s) + (dg^ -|- sin^ q dslg) (3.19) 

is given by 

z = ucosh.p r = 'usinhp (3.20) 

Setting u = e imposes a (dimensionful) holographic UV cut-off on the theory living on the 
boundary. Since the minimal value of the coordinate p is zero, it follows from (3.20) that 
the range of x is bounded hj z > e and hence the AdS^ cut-off p is related to the uniform 
UV cut-off e. 

At this point we do not have an analog of the map (3.20) for the bubbling solution that 
is valid for all values of the coordinates. We can construct the map for the asymptotic 
region dehned by the EG expansion for the AdS^ x boundary theory (2.21), however this 
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expansion breaks down once the FG coordinate u is not small. As discussed in [13] one can 
construct a map that is valid also near ^ = 0 by patching together the expansion near the 
AdS^ boundary and the FG boundary. We will not pursue this construction here since we 
focus all our calculations on the theory living on the AdS^ x S'^ boundary. However, since the 
metric is asymptotically AdS, we expect that one should obtain only a small modihcation 
to the identihcation of the UV cut-offs that, crucially, does not affect the logarithmically 
divergent term in the entanglement entropy. 

Physically, the interpretation of r; as a UV cut-off and the form of the subtracted en¬ 
tanglement entropy (3.18) is quite natural. Note that the dominant contributions to the 
entanglement entropy come from UV degrees of freedom located near the entangling surface. 
Since the Wilson surface (at hxed time) intersects the entangling surface at two points in our 
geometry (see hgure 1), the defect contribution to the entanglement entropy has essentially 
the same dimensionality as the entanglement entropy of a two-dimensional GFT. This is also 
reflected in the AdS^ slicing we employ and the fact that the minimal surface we hnd in the 
bulk (3.3) is familiar from AdS^/CFA 2 - Hence an argument along the lines of those given in 
[13] shows that the extra divergent contribution of the Wilson surface should be logarithmic, 
wherein the cut-off is associated with a minimal distance to the defect. 

4 Holographic stress tensor 

The goal of the present section is to calculate the one point function of the six-dimensional 
stress tensor holographically for an asymptotically AdSj x bubbling solution. Since the 
bubbling solutions are eleven-dimensional one has to utilize the machinery of KK holography 
that was developed in [39]. Note that a similar calculation was performed in [23] for the 
type HB bubbling solution dual to half-BPS Wilson loop defects and we will largely adopt 
their method to our case. 

The KK reduction of eleven-dimensional supergravity on produces a seven-dimensional 
supergravity (with negative cosmological constant) with inhnite towers of massive helds 
[40, 41]. These can be classihed by their seven-dimensional spin and representation of the 
relevant S'0(5) spherical harmonics (scalar and tensorial) on One difficulty is the mixing 
of modes coming from the eleven-dimensional supergravity as well as the fact that seven¬ 
dimensional fields can be related by eleven-dimensional diffeomorphisms, leading to nonlinear 
gauge symmetries. The resulting seven-dimensional action can be diagonalized and the 
masses of all the fields were determined in [40, 41, 42]. 

Via the AdS/GFT correspondence the seven-dimensional supergravity fields are dual 
to operators in the six-dimensional (2, 0)-theory. The precise dictionary can be found for 
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example in [43]. 

In [39] it was argued that in order to obtain a local seven-dimensional supergravity 
action without higher derivatives one needs in general to perform a KK reduction map that 
is nonlinear and relates the eleven- and seven-dimensionall helds schematically as 


= i’ll + +... 


(4.1) 


Here, /C is a differential operator and the ellipsis denotes higher order terms with three or 
more eleven-dimensional helds. This nonlinear mixing in general complicates holographic 
calculations (see e.g. [23]). However, a simple rule was derived in [39] to determine when 
and which nonlinear terms appear. For a supergravity held dual to a dimension A operator 
in the CFT, the only nonlinear terms that can appear are the ones for which the sum of 
the dimensions of their respective dual operators is less than or equal to A. When this rule 
is applied to the stress tensor, which has A = 6, it is clear that there can be no nonlinear 
mixing since the operators with lowest dimension have A = 4.^ 

The starting point for the calculation of the holographic stress tensor is to decompose 
the eleven-dimensional metric into a AdSy x part and a perturbation, denoted as and 
h respectively: 

= qmn dx^dx^ = + Hmn^ dx^dx^. (4.2) 


We use the FG coordinate chart (2.21) where we can identify 6 as the polar angle on The 
Wilson surface preserves an S'0(4) subgroup of the R-symmetry and so does the bubbling 
geometry. Therefore, performing the harmonic decomposition on the eleven-dimensional 
helds we obtain contributions only from spherical harmonics invariant under S'0(4). These 
depend only on the polar angle 6. The zero mode on of an eleven-dimensional held that 
only has nontrivial dependence on 6 can be expressed as 


(j){x) 


dO (j){x, 6) sin^ 6 
dO sin^ 9 


(4.3) 


The reduced seven-dimensional metric, which satishes the seven-dimensional linearized Ein¬ 
stein equation, is given by the following combination^ 


dsj = 


1 + -* 


<7‘“> + h 


flU 


dx^dx^ 


(4.4) 


^There is a ‘doubleton’ field dual to an operator of dimension A = 2, but such fields are free and decouple 
from the dynamics. 

^Our index conventions are: M, N,... are eleven-dimensional indices, ... are AdSj indices and a, 6,... 
are S'^ indices. 
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Here, is the AdSj vacuum metric whereas is the zero mode of the fluctuations in the 
seven dimensions. The held vf is the zero mode of the trace of the huctuations of the metric 
along the directions: 

7r{x,e) = (4.5) 

The factor of | in (4.4) comes from a Weyl rescaling to bring the KK reduced metric to the 
Einstein frame in seven dimensions. 

The calculation of the one point function of the stress tensor from the seven-dimensional 
metric (4.4) utilizes the standard method of holographic renormalization. In our case vf 
vanishes, and consequently the reduced metric is already in the EG form 

4/^2 

dsy = —— {dv? + gijdx^dx^) (4.6) 

where the large r limit corresponds to u —?• 0 and the metric gij can be expressed as a 
power series in u. The holographic stress tensor can then be calculated immediately using 
the formulae for d = 6 given in [44]. For the convenience of the reader and completeness 
we present the details of these calculations in appendix B. We have also checked that the 
counter-term approach developed in [45] gives the same result for the stress tensor. The 
hnal result for the change in the expectation value of the stress tensor in the presence of the 
Wilson surface on AdS^ x is 

A{Tij) dx^dx^ = (16 -h 3ml - Sms) {ds^^Ss “ dsls) (4.7) 

Note that the dependence of is as expected from a back-reacted supergravity solution. 
In addition, the expression depends only on the hrst two nontrivial moments m 2 , m 3 of the 
bubbling solution. The stress tensor contribution has a form that respects the so(2,2|]R) © 
so(4|M) of the Wilson surface and is traceless in line with the absence of a conformal anomaly 
on AdSs X S^. The use of this result is two-fold. First, any nontrivial holographic observable 
is useful to understand ( 2 , 0) theory better. Second, the stress tensor expectation value is an 
important ingredient in the calculation of the entanglement entropy using the replica trick, 
as we discuss briefly in 6 . 

5 Expectation value of the Wilson surface operator 

The expectation value for the Wilson surface operator can be obtained from the following 
formula: 


(Wr)=exp [-(/-/(o))] 


(5.1) 
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where / is the eleven-dimensional supergravity action evaluated on a general bubbling solu¬ 
tion and /(o) is the action evaluated on the AdS'y x vacuum. The action is given by 


/ = 


SnG 


( 11 ) 


N 


2 / ^ 


R - ^Fmnpq 


pMNPQj _ L J c AF AF 


+ / 


(5.2) 


The hnal term is the Gibbons-Hawking boundary term, which is necessary to make the 
variational principle well-dehned for spacetimes with boundary. The metric functions for 
the bubbling solution are given in (2.5) and the four-form held strength is given by 


F — {fl)^gim<^AdS3 A + {f2)^g2m^S^ A + {fsYgSmi^s^ A c"* (5.3) 


where e™ are the vielbeins on the Riemann surface S, ojx denotes the volume form for a 
unit-radius space X and the expressions for gj^n with / = 1,2, 3 can be found in appendix C. 


5.1 Action as a total derivative 

First we demonstrate the well-known fact that the on-shell action of eleven-dimensional 
supergravity is a total derivative. Using the Einstein equation, the Ricci scalar can be 
eliminated from the bulk term of the on-shell supergravity action (5.2): 

The equation of motion can be expressed as 

d* F + A F = 0 (5.5) 

The hrst term in the action can then be written as 

F A *F = dC A *F 

= d{C A*F) + C Ad* F 

= d{C A*F)-^C AF AF (5.6) 

Plugging this expression into (5.4) the C A F A F terms cancel and the on-shell value of the 
action is indeed a total derivative: 

/bulk =-^ [ die A *dC) (5.7) 
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The duals of the three contributions to the held strength are given by 


*F(i) 

*F(2) 

*F{3) 


/ 1 

^ 3^ f rdrbidO - debidr ) A A Ugs 


ft 

ftfl 


/I 


ftfl 


fi 


rdrb2d6 - deb2dr\ A UAdS^ A 


- rdrb^dO + -debsdr) A UAdSs A 


(5.8) 


Using the expression for the metric functions fi from (2.5), the ratios appearing in (5.8) can 
be expressed as 


/MsV hWf 

V /i ; 41U+(1-|G|2) 

ffifsV hWl 

\ f2 J 41U_(1-|G|2) 

[l^) — imr- 

Since the volume forms on the AdSs and the three-spheres are all closed, the bulk action 
(5.7) reduces to a total derivative over the two-dimensional Riemann surface S. In terms of 
the polar coordinates r, 0, which were introduced in section 2.1 by setting w = re*^, the bulk 
part of the action can be written as follows: 


7buik — 


1 

487rGSi^^ 

1 

487rG5i^^ 


d{C A *F) 

{drttr + dette) dr A d9 A UAdSs A A Ugs 


where ar and ae are given by 

IM! 

vt 

flA 


Ur = 


ae — 


^ rdr{b\) ^^^^^rdr{bl) -F ^^^■9^(63) 


2/: 


m 


2/1 r 


3 1 -pS -pS 1 -pS -pS 1 

' -Mbl) + ^iil-deibl) + ^-deibl) 


Vi 


Vi 


r 


(5.10) 


(5.11) 


We have shown that the on-shell action can be written as a boundary term: the bulk term 
reduces to a total derivative and the Gibbons-Hawking term is also boundary term. The 
boundary in the w plane has two pieces (see hgure 4). First one has the cut-off surface 
parametrized by rc{0,e), where e is the FG cut-off. Second one has the real line, which 
is parametrized in r,9 coordinates by r G [0, rc(0, vr, e)] and 6* = 0, vr, respectively. In the 
following we will evaluate the contribution from each piece in turn. 
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Figure 4: The two boundary components: the cut-off surface (red) and the boundary along 
a; = Re w (green). 


5.2 Gibbons-Hawking term 

The Gibbons-Hawking term in (5.2) is the boundary term that has to be added to the action 
in order to give a well-dehned gravitational variational principle. A ten-dimensional surface 
dehned by a constraint 


F (i“) = 0 


(5.12) 


can be parameterized by a set of coordinates a". The induced metric is dehned via 

dxM 


M N M 

lap = e„ ep guN, e„ 


da'^ 


and its determinant is given by 


\/l = ^ UJgs 


The normal vector (which we assume to be always space-like) is dehned via 

1 OF 

riM — 


dF^dll NP dx^ 
dx^ dx^^ 


(5.13) 


(5.14) 


(5.15) 


The extrinsic curvature and its trace are dehned via 


Ko^p = (Vm^tv) ej", K = (5.16) 

As discussed above the boundary has two components, which we now study in turn. 


5.2.1 Real line contribution 

The real line is dehned by F = y = 0, where w = x + iy, and the induced metric is simply 
obtained by dropping the Qyy component of the eleven-dimensional metric. The normal 
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vector is given by 


ft" = i (5,17) 

Using the form of the metric it is straightforward to determine the trace of the extrinsic 
curvature 




^dyfi ^dyf2 Sdyfs ^ 

2p/i 2pf2 2ph 2p2 


(5.18) 


The Gibbons-Hawking term along the real line can be determined from the expansion of 
the metric functions in the ?/ —)■ 0 limit given. This can be obtained from the formulae in 
appendix D. We do not present the details of this calculation but just present the hnal result 
which is 


J = \imVo\{S^)^Vo\{AdS 3 ) j dx [2Ay^ + 0{y^)] = 0 (5.19) 

which vanishes. This result was expected since at any point on the real line, one of the 
three-spheres shrinks to zero size and the space closes off. 


5.2.2 Large r contribution 

In this section we determine the contribution of the Gibbons-Hawking term from the large 
r cut-off surface dehned by the equation 


F{r,9) = r — rc{9,e) = 0 (5.20) 

for small e, where rc{9,e) is given by (2.26). Hence the surface extends along AdS^ and 
the two three-spheres and the induced metric in these directions is identical to the metric. 
We choose to use 9 to parametrize the cut-off surface and hnd the following nontrivial 
components of the vielbein e^\ 



ea 

^ 8 , 


el = de rc{9, e) 


(5.21) 


Using the formulae (5.15) and (5.16) we can calculate K and expand the result in a 
power series in e\ 




192 ^ 96m2 cos 9 


—3(16 — — dOms) -|- 5 cos 20(16 -|- 3m| -|- dOms) 

^ 8^2 

1 . o 

- cos 9 [— 8 m 2 -l- m 2 — 67712^3 -|- 10 m 4 

-I- cos 20 ( 2 dm 2 + ml — 10 m 2 m 3 -|- 70m4)] | + G 


(5.22) 
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After performing the integral over 6 we find that the Gibbons-Hawking term is determined 
completely by the contribution at large r, namely 


Igr = 


1 


L9 


^ d^^x^K= - ^Vo\{S^fVo\{AdSs) 


J SttGV 


5.3 Bulk supergravity action 


256 16 ^ , 2^ 


(5.23) 


The bulk part of the eleven-dimensional action integral (5.10) can be reduced to the following 
form: 


I = 




Vol(S''^)^ Vol(AdS' 3 ) dO dr {drttr + dgae) 


487rG 


Vol(S''^)"^ Vol(A(iS' 3 ) ® dr {hrar + flgag) 


N 


(5.24) 


«S(e) 


Here Vol(A(iS' 3 ) is the (regularized) volume of AdS^ since we consider the held theory on 
AdSs X and fii is the unit outward normal. If a boundary component is parametrized by 
r{T), 6 {T) then the unit outward normal vector is dehned as 


^ d9 dr 

Tir = — dr, ne = ——dr 
dr dr 


(5.25) 


There are three components to the boundary: 


dT,i : 0 G [0,7r], r = rc(0,e) 

8^2 : 9 = 0,r e [0,rc(0,£)] 
8^3 : 9 = n, r e [0, rc(7r, e)] 


(5.26) 


The contribution of c?Si corresponds to the cut-off surface, which we parametrize by 0 G [0, vr] 
and r = rc{9, e) as in section 5.2. The 9 dependence of r leads to an additional contribution 
to the normal vector 


I(cut) = 


VoliS^fVo\{AdS3) 




dr {hrar + hgae) 


9Si(£) 


1 f c)t 

-^ Vo\{Sy Vo\{AdS3) / d9(ar-^ae 

487rG!i^^ Jo V 99 


(5.27) 


r=rc(e,s) 


The integrations over 8 Tj 2 and 8 T ,3 combine to give the integration over the real line. Putting 
the two together and going back to the half-plane coordinates we hnd 


rrc(e= 0 ,e) ^ 

dr {hrar + hgag) = — dr - ag 

dT,2{£)+dT.3{e) Jo ^ 


rrcie=TT,e) ^ 

+ dr - ag 

0=0 Jo r 


(5.28) 


0 = 7 \ 
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Using ^de\ 0 =o = dy and ^de\e='K = —dy and the expressions for and ag given in (5.11), the 
contribution from the real line becomes 


1 |•r^{9=0,e) 

h-) = -^ Yo\{AdS^) / dx 


^/3 


2/3 2/3 


(5.29) 


5.3.1 Large r contribution 


Using the expansion of and ae at large r one hnds the integrand of (5.27) can be written 
as 


1 

Z9 




r=rc{0,e) 


5 sin 6 cos^ 6*(cos 20 — 5) (16 + 3m| — Sma) 

4e'^ 

sin 2^ 

H—[(2144m2 + 445m2 — 1760m2m3 — 860m4 

+ 20 cos 26 (64m2 + 15m3 — 807712^3 — 607774) 

+ cos 40 ( — 967772 — 257772 + 16077727773 + I4O7774)] 

+ O (e^) (5.30) 


After performing the integration over 0 the hnite terms above drop out and and we are left 
with 


I{cut) 




Vo\{S^fVo\{AdS3) 


487rGSi^^ 


-64 — 12777| + 327773 


+ 0(e^) 


(5.31) 


5.3.2 Real line contribution 

The expression from the real line can be obtained from expanding the integrand of (5.29) 
in the 7 / —)■ 0 limit. The behavior of the integrand in this limit depends on the interval x is 
located in. We dehne 


Jo = [-00, /i] U [^2, 6] U ■ ■ ■ U [/2n, +00] 
X+ = [ei, 6] U [^ 3 , ^4] U • • • U [e2n-l, 6n] 


(5.32) 


We present the details of the calculation in appendix D, with the hnal result given by 


lim 

J/-5-0 


f3f3 f3f3 f3f3 \ 

■^9y{bl) + ^-^dyibl) + ^-^i-dyibDj = 


16L^(2g^-rg3)(/>0 
( 35 ^- 92 ) ( 51 - 1 - 92 ) 

32L9(9f+39i92-93)(0O-l-23:) 
( 91 - 1 - 92 ) (9?+ 292 ) 


X eXn 


X e I4 


(5.33) 
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where the Qi and 0o are functions of x and can be found in (D.3) and (D.5), respectively. In 
this limit the integrand (5.33) is nonsingular for any hnite x on the real line. However, the 
integrand grows as x —)■ ±cx) and we can extract the divergent behavior coming from the 
region near the cut-off. We determine the divergent contributions in (D.14) and we write 
the result as 


I{x) — + 




Yo\{S^fYo\{AdS^) 




256 80 — 12m2 -l- 32m3 _ ^ / 2\ 

■7ir +--2 + ^ + o(£^) 


(5.34) 


The hnite term can in principle be determined by performing the x-integral for the full 
integration region and subtracting the divergent contributions given above, i.e. 


32(f?! + 3(71(72 — S's) (00 + 2x) f 16 { 2 gf + gs)(j)o 

(s? + 92)(s? + 292) b9?-92)(9? + 92) 

-(-f+ ^°-12,g + 32’"3 ^| (5.35) 

The last two terms remove the divergent contributions from the integral over Xq that is 
regulated for large positive x by x < Xc,+ (e) and for large negative x by x > Xc-{e). 

At this point we have not been able to hnd a closed expression for these integrals, but 
they can be evaluated numerically. If a relation to a matrix model calculation exists these 
integrals may be related to the resolvent. However, as no proposal for a matrix model exists 
at present we have not pursued the evaluation further and leave this for future work. 



5.4 Final result 


Combining all the contributions to the action, which can be found in (5.23), (5.31) and 
(5.34), we obtain for the on-shell action 


I = 


X9 


Vo\{S^YVo\{AdS3) 


StiG 

ATS 


( 11 ) 

N 


80 3 T ^ / 2\ 


Vol(AdS'3) ^ 16 ^ 


(5.36) 


The hnite contribution is given by (5.35). This can be evaluated exactly for the vacuum and 
the result is 


.^(0) = -64 (5.37) 

Thus, using (5.1) we can express our hnal result for the expectation value of the Wilson 
surface operator as 

log(iyr) = -^Vol(Ad^3) (.X + 64) (5.38) 

lyZTT 
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Note that the power divergences in e are independent of the details of the bubbling geometry 
and so cancel in the subtraction. However, the result is proportional to the inhnite volume 
of AdSs, which is regulated by the cut-off z = rj introduced in the entanglement entropy 
calculation in (3.13). 

As is the case for the stress tensor contribution, the expectation value of the Wilson 
surface operator constitutes a potentially useful holographic observable of the ( 2 , 0 ) theory. 
This holographic result may be compared to direct calculations in this theory as well as be 
applied to the replica calculation of the entanglement entropy. It would be very interesting to 
study localization and related methods to calculate the expectation value of Wilson surface 
operators in the future. 

6 Discussion 

In this paper we have calculated the holographic entanglement entropy for the six-dimensional 
(2, 0) theory in the presence of a Wilson surface. In addition we have calculated two other 
holographic observables: the one point function of the stress tensor in the presence of a Wil¬ 
son surface and the expectation value of the Wilson surface operator. The bubbling solution 
that provides the holographic dual of the ‘heavy’ Wilson surface is determined in terms of 
2 n real numbers that parametrize the partition of the real line into two types of segment 
(where G = 0 or G = - 1 -i). The quantities we calculated all depend on moments ruk whose 
definition in terms of is given in (2.18). In particular, the stress tensor only depends on 
m 2 and m 3 . On the other hand, the final term in the result for the entanglement entropy 
(3.18) cannot in general be expressed in terms of the moments mk- The finite part of the 
expectation value of the Wilson surface operator is expressed in terms of a real integral. 
While it is possible to evaluate this finite part for n = 1 and n = 2 in closed form, we have 
been unable to evaluate it in general. Some numerical experiments however indicate that 
this integral does not have a simple expression in terms of the moments m^ or the final term 
in the entanglement entropy (3.18). 

This ‘coloring’ of the real line is reminiscent of other bubbling solutions, such as those 
dual to Wilson loops [14] or the LLM solutions [46], where the data of the solution can be 
mapped to the data of a matrix model or free fermion phase space, respectively. One hint 
that BPS Wilson surfaces might be described by matrix models can be seen as follows. If one 
compactifies the six-dimensional (2, 0) theory on a circle one obtains five-dimensional SYM 
theory [47, 48] and a Wilson surface that wraps the circle becomes a Wilson line in this hve- 
dimensional gauge theory [49, 32]. Various quantities have been calculated for this system 
using localization and it would be very interesting to see if these results can be lifted to six 
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dimensions and compared to our calculations in the limit of very large ‘representations’, in 
analogy with [21]. 

It would also be interesting to see whether the calculation of the entanglement entropy 
in the presence of a Wilson loop in SU{N) W = 4 SYM due to Lewkowycz and Maldacena 
[17] generalizes to the Wilson surface in the (2, 0) theory. Recall that their calculation 
used the replica trick and involved the expectation value of the Wilson loop and the stress 
tensor in the presence of the Wilson loop on the space x H^. A generalization would 
most likely start from the expectation value the Wilson surface and the stress tensor in the 
presence of the Wilson surface on x H^. Since our holographic calculation gives these two 
quantities for the AdS^ x boundary, if it is possible to map the results to x then it 
should be possible to compare the Lewkowycz and Maldacena calculation to the holographic 
entanglement entropy calculation we have performed. One complication is that unlike a 
one-dimensional Wilson loop, the two-dimensional Wilson surface has a conformal anomaly 
[50, 51, 52] and it is not clear how to determine its contribution in our case. A simpler case 
in which to consider the anomaly might be the case of an abelian Wilson surface, as studied 
in [52]. 

Our calculations provide results that could be compared to any held theory or localization 
calculation in the (2,0) theory. What is missing at this point, compared to the analogous 
case of the half-BPS Wilson loop in A/" = 4 SYM, is a concrete dictionary between a matrix 
model and the supergravity parameters as well as a better understanding of the anomaly and 
the map of the stress tensor and the expectation value to the x boundary geometry. 
We leave these interesting questions for future work. 
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A Contributions to the entanglement entropy 


In this appendix we carefully discuss the contribution to the area integrals that are needed 
in section 3.2. 


A.l J\ 

First we consider Ji, given in (3.9). The radial integral can be performed directly, but it is 
useful to rewrite it in terms of Legendre polynomials using (2.15). We divide the integration 
range into two regions: 0 < r < |^j| and |.^j| < r < Tc {0,£)- For each region we choose the 
Legendre representation that converges, yielding 


2 n „ 

4=1 >^0 


f d9 sin 9 < 

rrc{9,e) r^(rcos6' —£ 

/ dr - 

^ Y. Pe (cos 9) ( 

/o [ 


£=0 ' 


“hd r^{rcose-^i) ^ ^ ^ 

dr -r^-;- > Pi (cos 0 ) 

e=o 


JO 1^*1 

Performing the two radial integrals directly we hnd 

..3 


(A.l) 


2 n „ 

J. = / 

i=l >^0 


d9 sin 6 * 


cos 0 Po y + (cos 9 Pi- Po) 6 y + (cos 9 P2- Pi) r 


+ (cos 9 Ps- Pa) log r - ^ 


(cos^P^+s - Pi+2) C 


i+o 


- rc (0, e } 


e=i 




16 


-Po 6 y + X] 


{cos9Pi_4- Pi_ 3 ) 


£=4 




£-4 




Orthogonality of the Legendre polynomials can be expressed via 

[ d9 sin 9Pi (cos 9)Pk (cos 9) = ^ 

Jo 2t + 1 

We use this to simplify the above expression dramatically: 

2 n ^ 

d9 sin^ 


(A.2) 


(A.3) 


J. = -4L>X;(-l)‘{ [ 

i=l '<^0 


3 2 

cos^Po ^ + (cos^Pi - Po)e* J + (cos^Pa - PO^-r. 


+ (cos 9 P 3 - Pa) 6 ^ log rc-Y^ 


(cos0P,+3-P£+2) 6 


£+3 


£=1 


2 .3 

H-£ 


(A.4) 


Note that the hnal term is a sum of contributions at r = |£j|. Substituting for the cut-off 
function rc{9, e) given in (2.26), we then expand in e up to and including O (£°) and perform 
the remaining integrals over 9. We hnd the hnal result 

64 —24 -|- 3m2 — 8177.3 


Ji = r 


3^4 


+ 


15 


+ 0(.^) 


(A.5) 
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A.2 t /2 

Next we calculate J 2 , which we reproduce from (3.11): 


Jo = —2L® / d6 sin^ 


/'rc(d,e) 


drr"^ 


X 


2 n + 2 '^(-1)*+^'_ - r cos6> _ 

i<j \/r^ - ‘2r^i cos d + - 2r^j cos 9 + C 


(A.6) 


We can split the integral into two terms coming from the sum in the last line of (A.6). The 
hrst term is simply 

4 


J 2 ,a = —^nL^ / sin 6* rc (6*, £)" 

O Jn 


= L« 


64 n 77,(40 + 3777n — 87773) ^ 

+ ^- ttA - — + O [e^] 


3^6 


18£2 


(A.7) 


The evaluation of the second term, denoted J 2 ,b = J 2 — J 2 ,a, is more involved than that 
of Ji. Our strategy is to divide up the radial integration range and replace the square root 
factors with the appropriate convergent series of Legendre polynomials in each interval. The 
fraction in the summand is symmetric under (7 -H- j) so we can choose |^j| < without 
loss of generality and write: 

_ p7T 

J 2 b = -4L® ^(-1)*+^' / de sin 6 

i<j -^0 




X 


r^(r^ - r cos e + Q) + 60) ^ f 




£ 00 


'hi I 


h^Jl £=o ''' '' k=0 

dr -r cose ^ ^ ^ 

e=o 


r\i 


3 1 


'0 


'1^*1 r^(r^ -r cos 9 + Q+ ^Q) ^ 


le^iio 


i =0 


k=0 
£ 00 


fc =0 


= K, + K2 + K3 


(A.8) 


where the Legendre polynomials are all functions of cos 6*, as before. First let us consider 
Ki. 

_ p7T 

ATi = -4L^ ^(-1)*+^' / de sin0 

i<j 


'0 




'•rc( 6 »,£) 


£,k=0 


'hd 


dr (r^ - r cos 0 r 


-£-k 


X 


(A.9) 












We must perform the radial integral first because its upper limit depends on 6 . This results 
in several sums over powers of r and a logarithm: 


Ki = 


_ °° p p c£ck „3—£—k 

-4L»5;(-1)‘+^ / ci9sin9.j ^ ^ 

i<j I e.,k=0 

£+/c7^3 


cos 


^(6 + 0 ) X] 


°° p p „2—£—k p p ^1—A; 


e,k=o 

t+k ^2 


2-i-k 






e,k=o 

t+k^l 


1-i-k 


/ 3 2 1 ^ 

^ p,p,efe^-cos0(e.+e,) PiPkiie.+iii, Y. logr 


£, fc =0 

V+fc =3 


£, fc =0 

£+ fc =2 


e,k=o 

£+fc=l 


/ 


rc( 0 ,£) 


hd 

(A.IO) 


We only require the entanglement entropy up to and including O (£°). Recall that the cut¬ 
off function rc{0,e) given (2.26) leads with 0 {e~‘^), and therefore only the logarithm and 
non-negative powers of r contribute to the upper limit. Specihcally, we can terminate the 
inhnite sums in the hrst and second lines at 3, 2 and 1, respectively. Integrating over 6 we 
hnd 


^upper ^ ^9 


64n n(40 -|- 3m| — Sms) 


3^6 


18£2 




where we have made use of the following results: 


(A.ll) 


^(_ 1 ).+; = - 0 )" = -4 (A. 12 ) 

i<j i<j 


Next let us consider the lower limit and perform the integral over 9. The terms with no 
explicit cos^ factor vanish unless £ = k hy orthogonality (A.3). To deal with the terms that 
do have an explicit cos 9 factor, let us dehne 

Xik ^ J sin 6* cos 6*P^ {cos 9)Pk {cos 9) = 2 ^ J ^ ^ ^ (A. 13) 

2{£-kf{l + £ + k) 

~ {£ + k){2 + £ + k){l + £-k)\{l-e + k)\ 

These terms are only non-zero when Xik is too, which occurs when \£ — k\ = 1. These two 
observations imply that the coefficient of the logarithm vanishes and that the conditions on 
the sums in the first two lines of (A.IO) have no effect for the lower limit. All that remains 
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IS 


i<j 


e,k=o 

(6 + ^j) 


2 e f\Q\ 

Oik 


2i + l 


3-e-k \c\l-t-k 

_+ M _f.f. 

i-k 1-i-k 


2-i-k 

Performing the sum over k and using the dehnition (A. 13) we hnd 


(A. 14) 


K'r’ = +4r’‘E(-i)‘*'E 5^44 

i<j £=0 


l?jl 


3-21 


\l-2l 




3-2i 1 - 2£ 


(2£ + 3)(l-2£) 


(A.15) 


The limits on the radial integrals in iP 2,3 are independent of 6 so we are free to reverse the 
order of integration. Let us begin with K 2 '. 


K2 = -4L^ Y1 


-l)*+i 




,k-l 


I I l^i| '^3 

pTT 

X de sin 6PePk (r^ - cos 9 + Q) + r^iQ) 


-4L»E 


(-l)*+i 


“IL-I 


k-l Si 


i<i 


fel 4|‘'bS' 4L2 «+i 


^ik (r^ + riiij) - r^{ii + ij)X, 


Ik 


(A.16) 


where again we have used (A.3) and (A. 13). Performing the sum over k then integrating 
over r we hnd 


^ l^jl ^2t+ll' 

i<j £=0 ^3 

^ ' A r, S*Si 


l-hlifl + li 

2£ + 3 ' * 


(A.17) 
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We can compute it's using the same method: 


i^j Isdisjl Jo 


X I d6 sin OPePk cos 6 {^i + ^j) + 
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- 4 i’E 


i<j 


, / ^^Y.Tifk 

jl ^'O gu=Q^l^j L 


le^iie 


2 £ + 1 


(r + r - r (^i + 


- 4 L»vW^f; 2 


1 


i<j 

3 + 2£ 


^=0 


a/ 


2^+ufe 




^+1 

2£ + 3 


K. + &)(^ + ^ 


J / J 




| 5 + 2 £ 


5-2£ 


(A.18) 


Now we combine the hnite contributions to J 2 ,b- The inhnite sums can be evaluated and 
the following remarkably simple result is obtained: 

4L9 


A'!”” + /n + A's = Yr ii (6 - &)' 

i<j 


(A.19) 


Recall that we have assumed |.^i| < \^j\. The ordering of the in (2.13) implies that sgn,^j 
evaluates to +1 and since the sum is ordered we can write the hnite contribution to J 2,6 as 

A/'S 

(a.20) 

i<j 

Note that this term cannot be expressed in terms of the moments mk- Thus, our hnal result 
for J 2 ,b is: 


J2,b — 




64n 


n(40 + 3m2 — Sma) 
l8£2 


64 

3^2 




(A.21) 


Summing (A.7) and (A.21) we hnd 


J2 



+0(A) 


(A.22) 


B Calculation of the holographic stress tensor 

In this appendix we present some details of the KK reduction calculation as well as the 
calculation of the stress tensor using holographic renormalization. As mentioned in section 4, 
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first one has to decompose the metric into the vacuum AdSj x S‘^ part and fluctuations, as 
in (4.2). In FG coordinates the vacuum metric is given by 


9^mn dx^dx^ = 




u u 


16 


— I du + I 1 H —;;— h — I ds + I 1 — ;; — 1“ “ 1 ds 


u u 


16 


’S3 


+d6‘^ + sin^ 6ds% 

03 


(B.l) 


The metric fluctuations in terms of the functions ai{u,d) appearing in (2.21) are 


hiuN dx^dx^ = L 




dsAdSs V ^ y ~ 16 ) 


+ (q !3 — 1 ) dO + ( Q !4 — sin 6 ) ds 


"S3 


(B.2) 


Using these expressions, we calculate the seven-dimensional reduced metric (4.4) and the 
outcome is 


dsj = 


4L^ 




u 


u 


du^ + ( 1 + + TF + ~ 8ms)w ds 


16 320 


’AdS3 


+ 1 


1 . 

- 1 -(16 -|- 3m: 

2 16 320 ^ 


8m^)u^ 


ds'gs 


(B.3) 


Notice that substituting the vacuum moments in (B.3) one can retrieve the AdSj entries 
in (B.l). This is because the trace shift does not contribute to the reduced metric, i.e. tt 
vanishes. Furthermore, a further FG map of (B.3) is not necessary since it is already in FG 
form: 

42,2 

dsl = —— idu^ -|- Qij dx^dx^) (B.4) 

where the six-dimensional metric gij is given by a power series in u\ 


9 = 9{o) + 9(2) + 9(A) + 5'(6) \ogu^ + ... (B.5) 

To compute the holographic stress tensor, we simply read off the asymptotic metric 
coefficients ( 7 ( 0 ), g( 2 ), g(A) and 77 ( 6 ) from (B.3) and substitute them into the d = 6 formula 
given in [44]. For completeness we present this fomula here: 

— g ^(7) (^9(6)ij — ^(6)ij + ^ Sij^ (B.6) 
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where the second and third terms are dehned via 




2(c/(2)c/(4))*i + {9{4)9i2))ij - {9i3))ij + ^9{2)ij - (tr t/(2))^) 

-tr^(2) (^9ii)ij - \{9l2))i}j -9m Qtr^(22)tr^(2) - ^(tr^(2))^ 
1 1 \ " 


6 


+ l^^^{9{2)9{i)) 


(B.7) 


Sij — - 2R^l-Cki + 4 ((5'(2)5'(4)) - (5'(4)5'(2)))j^- + - g^Q^ijV'^B) 

2 / 2 4 3 \ 

+ -gmB + 9{o)ij (-3 tr^f 2 ) - ^(trt/(2))^ +-trc/(2)trc/^2) j (B- 8 ) 

with the quantities Cij and B dehned by 

11 1 
Cij = gm - -{g‘(2))ij + ^5'(2)ji tr5((2) + -g(o)ijB 

B = tr gf^) - {tr g( 2 ) f (B.9) 

Note that the contraction of indices is performed with the inverse of g{o)ij- A general formula 

for the trace of the stress tensor follows from these dehnitions: 

3 (2L)^ / 1 1 1 1 \ 

(T\) = l^- — (trg(2)f + -trg^2)trg^2) " q^^9(2) + 3 tr(7(2)^/(4)J (B.IO) 

Evaluating these formulae we hnd 


/n^ \ ^ i 20 + Om^ 24m3 ^ 2 ^ 2 \ 

{Tij) dx dx^ = —^-—- [dsAds^ - dss^) 

oTTCj' 


(B.ll) 


Notice that the stress tensor is traceless, which rehects the fact there is no Weyl anomaly 
for AdS^ X S^. After observing 

1 voi(sy , 8^2 

and using the dehnitions (3.15), we then subtract oh the contribution from the vacuum to 
obtain our hnal result: 

A{Tij) dx'dx^ = (16 + 3m2 - Sms) {ds\^s^ - ds%) (B.13) 

C Four-form field strength 

In this appendix we present the formula for the four-form held strength 

F = {fl)'^gim^AdS3 A e™ -I- {f2)^g2m^S3 A 6^ + (/3)^5'3m A c”* (C-1) 
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where ojx denotes the volume form for a unit-radius space X. The gim are related to 
derivatives of potentials bj via 


Uifgiw = -d^bi/L^ = 2(j+ +jj 

{f2f92w = -d^b2/L^ = -2(i+ - i“) 

(C.2) 

Since the four-form field strength is related to the three-form potentials by F(/) = dC(^i), it 
follows from (C.2) that the potentials take the following form: 

C(i) = bi—dz Adt Adi 

( 7 ( 2 ) = &2 sin^ 6 i sin 62 dOi A 6*2 A dO^ 

(7(3) = 63 sin^ tfji sin'^2 d'ljji A d'02 A dtfj^ (C.3) 


Next we review the the expressions for the fields j in (C.2) found in [37]. The currents 
can be expressed in a compact way by dehning 

Jn, = [G (<7 - 3(7 + AGG^) d^G + G {G + G) d^G] (C.4) 

and are given by 


j+ = 2i ((G - Gf - AG^G) W-^ 

j- = 2 GJ^ (-2GG + 3G^ - G^ + AG^G^) W-^ 

f=3dh^ _2J + 

^ G(1 - GG) G(1 - GGf 


It is then straightforward to verify that the potentials are given by 


bi = 

b2 = 
bs = 


2{G + G)h 
2GG + i{G-G) 
2{G + G)h 


+ 2 h- 2 A> 


2GG - i{G 
{G + G)h 
'4(G(7- 1) ‘ 


G) 


<h 


2h + 2<h 


Here, h is the dual harmonic function to h and satisfies 

idyjh = dyjh 

With h = —iL^{w — w) as in (2.11), one obtains 


h = L^{w + w) 


(C.5) 


(C. 6 ) 


(C.7) 


(C. 8 ) 
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Also, $ is defined via 


Gd^h = (C.9) 

Using dwh = —iL^ and G given by (2.12) we solve (C.9) to find 

4 = L’5;(-l)A(tt.-«(©-{,) (C.IO) 

j 

Note that $ is real, hence the only thing that could be added is a constant, corresponding 
to an ambiguity in the dehnition of the 6/. 

D Calculation of the real line contribution to the on-shell action 

In this appendix we present details on the calculation of the contribution from the real line to 
the on-shell action. To do this we have to expand the metric factors and 6/ in a power series 
in y around y = Q. The important point is that the expansion of G, G differs in different 
intervals. Let us dehne 


Xq = [-(X), ^i] U [^ 2 , ^ 3 ] U ■ ■ ■ U [^ 2 n, + 00 ] 

X+ = [^ 1 , 6 ] U [ 6 , ^ 4 ] U ■ ■ ■ U [e 2 n-l, 6 n] (D-1) 


For the Taylor series expansion of G we have 

^ ^ f 0 gi{x)y + ig 2 {x)y‘^ g 3 {x)y^ + ... x elo 

1 i + gi{x)y + m{x)y‘^ + g3{x)y^ + ... x 

where 

«3(n = E(-it^‘iy4F 


(D.2) 


(D.3) 


For the calculation of bj we also need the Taylor series expansion of <F dehned in (C.IO): 


$ = {(j)o{x) + (t) 2 {x)y‘^ + (t)i{x)y‘^ + ...) 


(D.4) 
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where 


j 

j \ \ 

(D-5) 

^ I I 


Note that gi = 02 and = 204 which will be important in the expansion of the action. The 
combinations of metric fnnctions appearing in (5.9) can be expanded as follows: 


and 


and 



L^{9l-92)^ 

9 I +92 


+ 0{y^) 


L^{9f+92f 

2 ( 91 + 292 ) 


2/3 + 0(2/5) 


X eXq 
X eX^ 



+ :rel„ 



8L^ 1 

(9!-9i) 9® 



4L3(9?+2g2)^ 
9?+92 


2/3 + 0(2/3) 


X eXq 

X eX+ 


(D.6) 


(D.7) 


(D.8) 


The expansion of bj with / = 1,2,3 works the same way, bnt there are some cancellations 
due to the relations of the expansion coefficients for G and $ mentioned above. We hnd 


(^ - 2.#.„ + 4x + xeSj 

i0j^-2<#.o + 4a: + O(!/0 I £ 1+ 


(D.9) 


Since the subleading term is of order 2 /^ as 2 / —)■ 0 we hnd that dy{h\) is of order y. Also note 
that the metric factor (D. 6 ) is of order 2/3 as 2 / 0. Thus we hnd that the contribution to 

the action coming from hi vanishes at 2 / = 0 and hence does not contribute. 

The Taylor expansion of 62 is given by 


{ L0-j^ + 24+4i + O(!/0 I6l„ 

= < 

I L3 ((200 + 4a:) + {gl + 3 ^i ^2 - ^ 3 ) 1 /^ + 0{y^)) xeX+ 


(D.IO) 
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It is important to note that for x G X+ we find that dy{b 2 )^ will behave as a.s y ^ 0 and 
together with the behavior of the metric factor (D.7) produces a hnite contribution to the 
action. 

Similarly the Taylor expansion for 63 is given by 


f L3 (-00+(<?? +f)2/' +0(7/5)) 

63 = (D.ll) 


In a similar manner as for 62 we note that for x G Xq the dy{h‘iY term will be of order y^ 
which together with the behavior of the metric factor (D. 8 ) will produce a hnite contribution 
to the action at 7 / = 0. 

Summarizing we hnd that 

lim X G M 


2/1 


0 


32L'^(gf+3gig2-g3)((/’0+2a:) 

(9?+92)(fff+2g2) 

16L^ {2gl+g3)4>o 
(.9'i-g2){g'i+g2) 

0 


X G Xo 

X G X+ 
X G Xo 

X G X+ 


(D,12) 


The integration region Xo is cut off by the large Tc cutoff and includes the intervals [—Vc, .^ 1 ] 
and [^ 2 n, ?"c] that are responsible for Tc divergent terms. Using the large |x| expansion one can 
show using the Taylor series expansions of (D.12) for large arguments that the contribution 
from the integral is given from the large integration limits Xc,+ and Xc- by 


dx 




r-a 


m 

2/1 

p3 x3 


dy(h\)\y=Q = X® [—Ibx^^ + 12 m 2 X^^^ + (16 — + 16 m 3 )xc,+] + hnite 


fif'd 

I ^y(^Dly=o ~ [~16|xc,_|5 — 127712|xc_P + (16 — 9ml + 167773 )|xc,_|] + hnite 

'xc- 2/3 


(D.13) 


Using the fact that Xc,+ = rc(0, e) and Xc- = rc(7r, e) together with the relation of the radial 
cut-oh to the FG cut-oh parameter £ given in (2.26), one can extract the contributions of 
the x-integral that are divergent with respect to the cut-oh e: as follows: 


rdiv 


X9 


487rGSi^^ 


VoliS-^yVoliAdSs) 


256 , 80-12m| + 32m3 


(D.14) 
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